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An exact-diagonalization technique on small clusters is used to calculate the dynamical density 
correlation functions of the dimerized t-J chain and coupled anisotropic t-J ladders (trellis 
lattice) at quarter filling, i.e., the systems regarded as a network of pairs (dimers or rungs) of 
sites coupled weakly via the hopping and exchange interactions. We thereby demonstrate that 
the intersite Coulomb repulsions between the pairs induce a low-energy collective mode in the 
charge excitations of the systems where the internal charge degrees of freedom of the pairs play 
an essential role. Implications to the electronic states of NaV^Os, i.e., fluctuations of the valence 
state of V ions and phase transition as a charge ordering, are discussed. 
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A compound ev'-NaV^Os discovered as an inorganic 
'spin-Peierls' systernEP has attracted much attention be- 
cause of its intriguing charge degrees of freedom related 
to fluctuations of the valence state of V ions. It has been 
reported, e.g., that the ratio of the size of the spin gap at 
T = K to the transition temperature is largely deviated 
from the BCS. mean-field ratio as. 2A/fceTc = 6.5 with 
A = 9.8 mevBS and T c = 34 rijf that a large amount 
of entropy is released ai the phase transition in contrast 
to the case in CuGeOaJSllfand that anomalous behavior 
of the dielectric constantoLp along the a-axis occurs at 
T c while no anomaly is observed along the b-axisfP Be- 
low T c , a huge increase in the thermal conductivity has 
been reported to occur,® and with a direct evidence, 
the charge-ordering scenario for the mechanism of the 
phase transition has been put forward from NMR ex- 
periment.© Some possible ordering patterns ^ -charges 
at T < T c have been suggested from theories.E3E3li3'E3 
It has also been pointed out that r the-.<~l eV features of 
the observed optical conductivity!!^ E-§ cannot be inter- 
preted without assuming the .-system to be charge dis- 
proportionated even above T c Ji5li3) and to be reconciled 
with the picture of the uniform valence^f_y_ir>ns which 
X-ray structural analyses have shown JUtES V3 we have 
argued that the rapkLcharge oscillation should exist in 
the system above T C S3 

There are also a number of experiments detecting fluc- 
tuations above T c . A recent line-shape analysis in the 
X-ray diffuse-scattering experimentEU suggests the pre- 
transitional structural fluctuations to appear, where the 
correlation lengths start to increase around ~50 K, ^70 
K, and above 90 K for the c-, a-, and fo-axes, respec- 
tively, and diverge as a power low by lowering tem- 
perature down to T c . There has been an independent 
observation of the presence of critical scattering from 
soft phonons above T c £3 Such a pretransitional effect 
has also been noticed in the ultrasonic experiment ,B3) 
where the anomalous behavior of sound velocity starts 
to develop at ^80 K. Anomalous increase in the oscil- 



lator strength of the low-energy continuum observed^ 
T ~ 150 K in the optical conductivity experimented 
might also be related to this effect. All of these exper- 
imental data seem to suggest that to clarify the mecha- 
nism of the phase transition of this compound it should 
be of primary importance to study the nature of the fluc- 
tuations and ordering of electronic charges; we therefore 
consider this issue in the present paper. 

The high-energy electronic states of the V-0 network 
of the (V2O3) 3 " 1 " layer may be described by the d-p model 
with strongly-correlated d^y-orbitals of V ions bridged 
by the p x - and pj,-orbitals of O ions. Our analysisE3 
has shown that the low-energy states of the d-p model 
may be mapped to the t-J ladder model with strong 
anisotropy between legs and rungs. There exist interac- 
tions between ladders, and thus the generic model for the 
low-energy electronic states is the coupled anisotropic 
ladders described by the trellis-lattice t-J model at quar- 
ter fillingO This system is a Mott insulator because of 
the strong anisotropy of the laddersO with increasing 
anisotropy the rung turns into an effective, jingle site 
at half filling with a Hubbard repulsionSKSc^P How- 
ever, because the anisotropy is not extremely strong, the 
system cannot be regarded as the one-dimensional (ID) 
Mott insulator, as is evident in the significant deviations 
of the calculated single-particle spectral functions for the 
ladders with^prxesponding anisotropy from those for the 
single chain JH3'l3 although the spin degrees of freedom 
behave as a ID Heisenberg chain ssath. Xhc effective anti- 
ferromagnetic exchange coupling.ISEZrEf The rungs (or 
dimers) with their internal hopping and exchange inter- 
actions are thus the essential ingredient for the basic elec- 
tronic states of the present system. A rung has two sites 
with one electron and is uniaxial in structure, and be- 
haves as a 'repulsive-?/ center' with the effective Hub- 
bard interaction U e g — It — J (for a t-J dimer)jH3) which 
however has an internal charge degrees of freedom; i.e., 
the one-electron wave function is a linear combination 
of the two tight-binding basis functions located on the 
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left and right sites of the rung. The rungs are aligned in 
parallel on the two-dimensional (2D) plane to form the 
trellis lattice, and this configuration produces a unique 
interplay between charge and spin degrees of freedom, 
which we believe is the source of the observed anoma- 
lous behaviors of NaX^Os. Here the coupling between 
the rungs, especially the intcrsitc Coulomb repulsions, 
plays an important role in the charge-ordering .scenario 
for the mechanism of the phase transitionO'EJ'H 

In this letter, we address the following question: sup- 
pose an assembly of such dimers and assume that the 
interaction between dimers is essentially the long-range 
Coulomb repulsion, then what is the charge dynamics 
of the assembly? We take two lattice models shown in 
Fig. 1: the 2D trellis lattice consisting of anisotropic 
ladders as a model for the (V203) 3+ layer of NaV20s 
and the ID lattice with dimerization as a model for the 
rungs coupled in the direction perpendicular to the lad- 
ders of the trellis lattice. We will thereby show that the 
intersite Coulomb repulsions between the dimers induce 
a low-energy collective mode in the charge excitations of 
the systems which is due to the internal charge degrees 
of freedom of the dimers. We will then consider its rel- 
evance to the electronic states of NaV205 and present 
some consequences on the low-energy charge dynamics 
of this compound. i_. 

The Hamiltonian of our model systems is given bytr 



cell) ring for the ID dimerized lattice model. 
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where c\ a — c\ a (l — ni- a ) is the constrained electron- 
creation operator at site i and spin a (=T, I), = 
n^+riii is the electron-number operator, Si is the spin-i 
operator, and (ij) represents a pair of neighboring sites i 
and j. We restrict ourselves to the case of quarter filling. 
We define the hopping and exchange parameters Uj and 
Jij as t± and J± for the rungs and t\\ and Jn for the legs, 
respectively, and also as t xy and J xy for the zigzag-chain 
bonds where the intersite Coulomb repulsions V xy is also 
taken into account. For the ID chain model, we define 
the parameters for alternating bonds as tj_ and J± for 
stronger bonds and t xy and J xy for weaker bonds. We 
include V xy only in the latter bonds. We use values of 
the ladder parameters obtained in ref.:li3 t± — 0.298, 
t\\ = 0.140, J± = 0.0487, and J\\ = 0.0293 in units of 
eV. The value t xy = 0.05 eV is assumed and the relation 
Jxy — 4:t X y/Ud is used for simplicity to obtain the value 
of J xy . The value of V xy is then varied for simulating 
various situations. The repulsion Vji on the leg bonds is 
also taken into account for some discussions. We employ 
a Lanczos exact-diagonalization technique on small clus- 
ters to calculate the ground state and excitation spectra 
of the systems. We take a cluster of coupled four 2x2 
ladders (or 2x2 unit cells) with periodic boundary con- 
dition for the trellis-lattice model and a 16-site (or 8 unit 



Fig. 1. Schematic representation of the lattice structures of (a) 
the trellis lattice t-J model and (b) the ID dimerized t-J model. 
The rungs in (a) and dimers in (b), which have the strongest 
hopping and exchange interactions, are indicated by bold lines. 



A direct way for examining the charge dynamics of 
the systems may be to calculate the dynamical density 
correlation function defined by 

C(q,u) = KV^M^o}| 2 S(w - (E v - E )), (0.2) 

where p q is the Fourier transform of the electron-number 
operator, and E v and \tp v ) are respectively the z^-th 
eigenenergy and eigenstate of the system (with the 
ground state denoted by v = 0). We have examined 
the individual excitations of electrons by the optical con- 
ductivity calculations,!^ whereas by the present C(q, uj) 
calculations one can find possible low-energy collective 
excitations as well. 

Now let us first examine the calculated results for the 
ID dimerized t-J model. The results are shown in Fig. 2, 
where the parameter values used are in the region of 
strong dimerizatioLL.2:_L Aa:;/ = 5.96 and thus the charge 
gap clearly opens.E3 We immediately find here that at 
V xy = the spectra are the ones consistent with particle- 
hole transitions across the charge gap but with increas- 
ing the value of V xy there appears the low-energy peak at 
q = which is smoothly connected to the higher-energy 
peaks at larger momenta with a well-defined dispersion 
ojq, the width of which scales with V xy . We thus find a 
collective-mode-like behavior in the spectra. Although 
not related to the present issue, we find for very large 
values of V xy ( > 1 eV) that there appears a set of low- 
energy peaks at w<0.2 eV, the position and weight of 
which are almost momentum-independent; this feature 
becomes dominant for V xy > 2 eV and is simply a conse- 
quence of the stabilization of the state with two electrons 
in a dimer, as is naturally expected by intuition. 

Similar spectral features are also found in the trellis lat- 
tice model (see Fig. 3) although the results are only for 
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Fig. 2. Dynamical density correlation function C(q, u>) calculated 
for the ID dimerized t-J model with the parameter values of t j_ , 
txy, J± t and J xy listed in the text. All the <5-functions in Eq. (2) 
are Lorentzian broadened with the width of 0.03 eV. 



the momenta q = (0,0), (it, 0), (0, tt), and where 
the spectra for the latter two momenta are equivalent due 
to symmetry. We find that the features reflecting indi- 
vidual excitations at small V xy values gradually change 
into the features resembling a collective-mode excitation 
by increasing V xy values, i.e., the low-energy peak at 
q = (0, 0) which 'connects' with the higher-energy peak 
at q = (tt, 0) and (0, 71"), where the width of the dispersion 
ujq again scales with V xy . If we also include the repulsion 
V\\ in the legs of the ladders, we find that the lowest- 
energy peak shifts to the nonzero momenta q = (0, 7r) 
and (tt, tt) as shown in Fig. 4, the behavior being nat- 
urally expected from simple consideration of the repul- 
sion between two electrons. Thus, either 'ferroelectric' or 
'antiferroelectric' fluctuation with corresponding disper- 
sion ujq is realized depending on the relative strength 
of the intersite repulsions; other types of fluctuations 
would also be realized if we include some longer-range 
repulsions. These observations are quite natural because 
the parameters we are working with are in the region 
where the long-range charge order appears in the infi- 
nite systems, as the calculated values of the equal-time 
charge correlations (niTij) directly indicate, (where there 
should be a critical repulsive strength V^JIiJ and thus we 
should observe the features corresponding to the Gold- 
stone mode at V > V c , although in finite-size systems 
no real symmetry-breaking occurs and thus there is no 
exactly zero-energy excitation. 



Fig. 4. As in Fig. 3 but for the case with nonzero values of the 
repulsion Vy . 



Fig. 3. Dynamical density correlation function C(q, w) calculated 
for the trellis-lattice t-J model. The momentum q = (q±,q\\) is 
shown in each panel. 



A speculation to what happens at T > in the ther- 
modynamic limit seems to be straightforward from the 
results given above. With increasing temperature from 
K the thermal fluctuations in the long-range charge- 
ordered state develop and at T = T c there occurs a 
phase transition. The fluctuations of charges would be 
remnants of the quantum fluctuations obtained for small 
clusters at T — K where the charges oscillate coher- 
ently between left and right sites of the rungs using their 
internal degrees of freedom. The phase transition would 
be of a classical order-disorder type as in the theory of 
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refB where the result of a 2D Ising model was used, 
and would thus be of the second order. Here we should 
note that the effective local-field and thus local order pa- 
rameter for the charge ordering persists even well above 
T c as we haye demonstrated via the optical conductivity 
calculations^' EJf where we used the term 'charge dispro- 
portionation'. In other words, we have the local order 
parameter as the spatial left/right shift of the electron 
in the rung (which is nonvanishing even at high tempera- 
tures) and its thermal (and quantum) fluctuations char- 
acterize the electronic state of the system. This is thus 
quite unlike, e.g., the usual charge-density- wave transi- 
tion where the normal fermi-liquid is recovered immedi- 
ately above T c and the location of the fermi momentum 
k-p plays an essential role. In the present system, the 
large charge gap exists, and there is no metallic screen- 
ing of long-range Coulomb interactions. At tempera- 
tures far above T c , say at room temperature, the local 
order parameters would be in a state of rapid and ran- 
dom oscillations without coherence, but with lowering 
temperature, the spatial correlation develops, and the 
correlation length increases and diverges as the temper- 
ature approaches T c where the long range order comes 
out. Fmnx^.Ynprimental data which we have summarized 
above ,E3 : E2l'ESiy one may expect that the temperature 
at which the coherent fluctuations start to develop would 
be in a range of T ~ 50 — 150 K. The character of the 
fluctuations, i.e., ferroelectric, antiferroelectric, or what- 
ever, would have a consequence on the charge-ordering 
pattern at T < T c , and needs to be clarified further from 
both experiment and theory. We should note that the 
lattice degrees of freedom may well play a role in the 
fluctuations of charges. A recent study of the local Hol- 
stein and/or Peierls couplings in the ID chain and ladder 
modelscS' seems quite suggestive, although nonadiabatic 
treatment will be essential for treating the slow order- 
parameter fluctuations with a time scale compatible with 
the lattice fluctuations. As for the spin degrees of free- 
dom, it is on one hand an open question whether the uni- 
form susceptibility in_such a charge state fits well to the 
BonnepFisher curveES with significant deviations below 
250 K,c3) but on the other hand the charge-ordering pat- 
tern would natucally explain the origin of the spin gap 
at T < T c .yiiiy) Actual determination of the pattern 
would however be an intricate problem because not only 
the long-range nature of the Coulomb interaction but 
also the lattice degrees of freedom seem to play an im- 
portant mleHS* Suggestions from experiment are most 
desirable Finally, let us remark on the submillimeter- 
wave ESR experiment by which a direct observation of 
the spin gap has recently been made successfully!^ We 
here would like to suggest that if the mechanism of the 
spin-gap opening is due. Id _bcal spin-singlet formation 
by the charge orderingBt!3E3 and if the slow fluctua- 
tions leading to the charge-ordering instability present 
well above T c , the fast time-scale measurement like this 
ESR experiment could possibly detect the 'spin gap' even 
above T c which is continuous to the spin gap below T c . 
So far, neutron scattering experiments have not found 
such signals. BtP 

In summary, we have calculated the dynamical density 



correlation function of the ID dimerized t-J model and 
trellis lattice t-J model at quarter filling using an exact- 
diagonalization technique on small clusters, whereby we 
demonstrate that the intersite Coulomb repulsions be- 
tween the dimers (or rungs) induce a low-energy collec- 
tive mode in the charge excitations of the systems where 
the internal charge degrees of freedom of the dimers (or 
rungs) play an essential role. We have argued that the 
electronic states and charge ordering of NaA^Os should 
be understood in terms of the fluctuations of the local 
order parameters representing the electron locations in 
the rungs, and have discussed the experimental conse- 
quences of the fluctuations. 
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